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"The  report  appears  to  solve  the  problem  of 
predicting  in  advance  of  experiment,  the  velo¬ 
city  of  detonation  of  an  explosive  substance, 
if  its  huat  of  combustion  and  density  are  known  - 
with  an  accuracy  equal  to  the  average  experi¬ 
mental  determination.  This  result  seems  to  me 
to  be  of  the  greatest  scientific  importance, 
and  it  is  hoped  that  the  substance  of  this  re¬ 
port  may  ultimately  be  published  vjhere  it  will 
be  available  to  scientists,  and  where  full  credit 
may  be  given  to  Professors  Kistiakowsky  and 
Wilson. " 


(5)  Five  copies  forwarded  to  Dr.  Irvin 
Stewart,  Secretary  of  the  National  Defense  Research 
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G.  B.  Kistiakowsky  and  E.  Bright  Wilson  (tr'55)  Offi¬ 

cial  Investigators,  Division  B,  Serial  No.  6. 
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by  this  method  of  prediction,  using  heats  of  combus¬ 
tion  and  additional  observations  to  check  theory 
being  obtained  as  part  of  NDRC 
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The  hydrodyne.mlc  theory  of  Part  I  should  permit  the  calcula¬ 
tion  of  the  detonation  velocity  of  an  explosive  if  the  following 
information  were  available:  composition,  heat  of  formation, 
density,  and  equation  of  state  of  the  products.  Unfortunately 
no  completely  satisfactory  a  priori  equation  of  state  data  are 
available  for  these  high  temperatures  and  pressures.  However, 
an  approximate  form  of  the  equation  of  state  (Becker)  is 
available  containing  one  unknown  constant  (which  is  a  function  of 
the  composition  of  the  products,  and  of  the  temperature,  but  is 
independent  of  pressure).  This  constant  is  of  the  nature  of  co¬ 
volume.  If  this  eovolume  could  be  calculated  for  a  given 
composition  of  the  products,  the  detonation  velocity  of  a  new 
explosive  could  then  be  obtained  using  only  known  or  readily 
obtainable  data,  in  most  cases. 

In  this  investigation  the  covoluoes  were  obtained  by  working 
backward  from  the  measured  detonation  velocities  at  varying 
densities  of  loading  for  a  considerable  number  of  explosives. 
These  results  (given  in  Table  4,  where  k^  is  the  covolume  in 
question)  show  first  that  kij  does  not  vary  much  with  density 
of  loading,  which  confirms  the  assumption  that  the  form  of  the 
equation  of  state  used  is  approximately  correct.  Secondly,  the 
covolumes  obtained  in  this  way  can  be  expressed  as  a  simple 
function  of  the  compcsiticn,  namely  as  an  additive  function  of 


the  covolumes  of  the  constituent  molecules  of  the  product.  The 
values  assigned  to  the  covolumes  of  the  Individual  species  vere 
adjusted  empirically  so  that  the  covolumes  of  the  mixtures 
agreed  vlth  those  found  as  above  and  listed  In  Table  6.  Actually 
for  most  of  the  organic  explosives  only  four  different  covolumes 
occur  so  that  this  whole  procedure  amounts  essentially  to  fitting 
the  observed  detonation  velocities  to  a  formula  which  Involves 
the  heat  of  formation,  the  density,  the  chemical  composition  and 
four  empirical  constants.  Furthermore  these  constants  agree 
suirprlsingly  well  with  the  van  der  VTaals  covolumes,  when  reduced 
to  the  same  temperature. 

Having  succeeded  In  reducing  a  considerable  body  of  data  In 
this  way,  It  becomes  possible  to  reverse  the  process.  The 
empirical  specles-covolumes  obtained  above  may  be  used  to 
calculate  detonation  velocities  for  new  explosives  as  yet  un¬ 
measured.  Furthermore,  this  calculation  can  be  carried  out 
approximately  for  substances  which  have  not  been  prepared  at 
all,  because  the  necessary  heat  of  formation  and  density  can 
be  estimated  from  the  properties  of  related  compounds  emd  from 
various  additivity  rules,  etc. 

The  results  described  above  are  tabulated  In  Tables  4,  5, 
and  6.  The  text  of  this  report  discusses  the  method  of  calcula¬ 
tion  and  the  nature  and  effect  of  the  approximations  and 
assumptions  used.  Appendix  II  gives  in  concise  practical  form 
the  procedure  for  carrying  out  the  numerical  calculations  in 
any  given  case. 
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Introduction 

In  a  preliminary  report  under  the  same  title  (hereafter 
referred  tc  as  Part  I),  the  authors  discussed  the  general  theory 
of  detonation  and  shock  waves  as  developed  by  Hugoniot;,  Raleigh, 
Jouguet,  Becker  and  others.  The  present  report  gives  an  account 
of  the  theoretical  work  done  by  Dr.  R.  S,  Halford  in  cooperation 
with  us  in  an  attempt  to  develop  theoretically  rigorous  and  yet 
practically  useful  equations  from  which  the  detonation  rate  (and 
Incidentally  the  pressure  and  temperature  in  the  detonation  wave) 
could  be  calculated  for  any  explosive.  The  only  data  which  the 
theory  requires  from  experiment  are  the  composition,  heat  of 
formation,  and  density  of  the  explosive. 

This  goal  has  not  been  quite  reached  as  yet  but  none-the- 
less  the  present  form  of  the  detailed  theory  enables  one  to 
calculate  the  detonation  rate  of  most  explosives  (i.e.  all 
organic  materials  and  some  others)  with  very  few  ^  hoc  assump¬ 
tions  and  with  an  accuracy  that  compares  favorably  with  the 
accuracy  of  experimental  measurements  of  this  quantity.  The 
prediction  of  the  detonation  rates  of  new  compounds  becomes 
therefore  possible  and  if  the  accuracy  desired  is  not  too  great, 
this  can  be  done  without  any  experimental  work  at  all,  because 
the  heats  of  formation  and  the  densities  can  be  estimated  from 
comparison  with  known  analogous  compounds.  Each  calculation 
involves  not  more  than  20-30  minutes  of  arithmetical  work  and 
hence  ma.y  be  strongly  recommended  to  all  who  intend  to  synthesize 
new  compounds. 

At  the  end  of  the  report  we  have  appended  a  brief  statement 
of  the  procedure  to  be  followed  in  calculating  the  detonation 


velocltjr  of  a  nev  compound.  This  is  intended  for  those  vho  do 
not  want  to  read  through  the  rather  extensive  mathematical 
operations  of  the  main  body  of  the  report. 

Certain  limitations  of  the  theory  must  be  stressed  now. 

The  theory  gives  only  the  greatest  possible  detonation  velocity 
under  given  conditions  of  loading  and  hence  presumably  the  max¬ 
imum  brlsance.  Ifnenever  the  initiation  of  detonation  has  been 
Inadequate  and  detona.tlon  is  chemically  Incomplete,  less  vigorous 
action  must  be  expected,  No  information  as  to  the  ease  of 
initiation,  that  is  the  sensitivity  of  the  explosive,  is  to  be 
gained  from  the  theory.  Finally,  the  theory  breaks  down  when 
one  of  the  products  is  partially  gaseous  and  partially  solid  or 
liquid  under  the  conditions  existing  in  the  detonation  wave. 

Only  when  each  product  is  present  substantially  quantitatively 
in  one  phase  will  the  theory  give  a  correct  result. 

\  1.  General  theory  applied  to  solid  explosives 

In  Part  I  the  following  general  equations  were  derived; 


where  D  and  w  denote  the  detonation  velocity  and  the  material 
velocity  in  the  wave  respectively,  the  subscripts  1  and  2 
denote  quantities  relating  to  the  state  before  the  passage  of 
the  wave  and  to  the  state  existing  immediately  behind  the  crest 
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of  the  -wave,  P  is  the  pressure,  p  the  density  and  E  the  energy 
(thermal  and  cheraical)  per  gram  of  material.  The  follovlng  equa¬ 
tion  defining  the  detonation  velocity  was  also  obtained  in  Part  I 
■(see  Eq.  1-5-6  and  1-5-7): 

w 

c 
O 


The  derivative  is  to  be  taken  adlabatlcally  and  reversibly,  v^ 
and  ^2  '^^^®  volumes  o.ccupied  by  M  grams  (for  instance  1  mole) 

of  the  explosive  before  the  passage  of  the  wave  and  by  the  detona¬ 
tion  products  derived  from  M  grams  in  the  wave.  It  is  assumed 
that  the  chemical  reactions  have  been  completed  by  the  time  the 
crest  of  the  wave  reached  a  given  spot  in  the  explosive. 

In  solid  explosives  the  Initial  pressure  is  of  the  order  of 
one  atmosphere  whereas  ?2  Atm.  and  greater;  vg  on  the 

other  hand  is  only  slightly  smaller  than  vi  and  therefore  P^  and 
P^v^  will  be  neglected  in  the  following  compared  with  ?2  ^2^2* 

The  resultant  error  is  wholly  negligible  but  the  simplification 
is  great  because  the  equa.tion  of  state  of  the  explosive  before 
detonation  is  eliminated  from  consideration.  For  the  gaseous 
products  of  detonation  the  following  rather  general  equation  of 
state  will  be  introduced  now: 

flK  \2 


(5) 


where  fd4,7jis  an  unspecified  function  of  the  volume  of  M  grams 
and  of  temperature  and  n  is  the  number  of  moles  of  gaseous_ 
products  formed  from  M  grams  of  the  explosive.  The  function,  f. 


is  also  a  function  of  the  mole  number  n.  In  general  n  will  change 
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with  T  and  P  "because  of  shifts  In  the  various  equilibria,  thus 
causing  f  to  depend  Impllcltljr  on  T  and  v  through  n  as  well  as 
explicitly  on  T  and  v.  In  what  follows  this  implicit  dependence 
will  be  ignored,  loe.  n  treated  as  a  constant.  The  more  general 
treatment  introduces  but  slight  gain  In  accuracy  as  will  be  seen 
and  produces  great  complication. 

Substituting  the  equation  of  state  (5)  into  the  well-known 
thermodynamic  relation: 


one  obtains  for  the  velocity  of  detonation  the  new  expression: 

where  f’(v)  stands  for  the  partial:  t  —  ^jhe  quantities 

w,  Pg  and  D  are  now  eliminated  from  equations  (2)  and  (3)  giving 


the  following  expressions: 


Cp  _ 


Cv 


(8) 


_  nP  Tk 


n(E,-E,)-- 


(9) 


Here  C  stands  for  the  true  heat  capacity  of  the  gaseous  mixtures 
under  the  conditions  of  the  detonation  wave.  The  energy  change 
in  equation  (9)  ma,y  be  replaced  by  the  following  expression 
derived  from  the  First  Law  of  Thermodynamics: 


+  W(T2-T,)Cv.- 


(10) 
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vhere  is  the  heat  absorbed  (Q  is  positive  vhen  heat  is  absorbed) 
in  the  reaction  of  M  grams  of  reactant  to  give  the  assumed  conpcsi' 
tion  of  products  and  is  mea.sured  in  a  constant  volume  calorimeter 
at  the  temperature  of  the  explosive  (T^,  normally  the  room  tempera¬ 
ture).  Cyi  is  the  mean  heat  capacity  (from  to  T^)  of  M  grams  of 
the  product  mixture  divided  by  the  total  number  of  moles,  N  (in 
the  ideal  gas  state  or  as  solids  or  liquids  at  low  pressures). 

The  coefficient  (dE/dv)fp^  must  be  evaluated  from  the  equation  of 
state  of  the  product  mixture  using  conveniently  the  relation: 


Until  now  the  calculation  has  been  carried  out  rigorously, 
but  the  results  are  too  cumbersome  for  practical  use. 

Extensive  studies  have  revealed  that  a  very  convenient  method 
for  the  solution  of  the  problem  is  to  attack  it  in  two  steps.  The 
first  neglects  gas  imperfections  and  introduces  certain  other 
simplifications,  while  the  second  either  eliminates  these  approxlms 
tlons  or  demonstrates  their  unimportance  for  the  final  result. 


2, _ Calculation  of  ’’ideal"  detonation  velocity 

Replacement  of  the  equation  of  state  (5)  by  P^v  =  nj_RTj_ 
where  n^^  is  the  number  of  moles  of  gaseous  products  leads 
immediately  to  the  following  simplification  of  equations  given 
previously:  (the  subscript  i  refers  to  the  "ideal"  state) 


V?  __  ^p,-  Y 

-- V I  {Tz'^  ~1T )  ■“ 


rnR  k: 
■  2n 


n;T?  lay 

rr  / 


(12) 

(13) 


(14) 
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Prom  these  equations  the  detonation  velocity  Is  readily 
calculated  when  the  composition  of  the  product  mixture  is  known. 
The  heats  of  reaction  are  calculated  hy  the  First  Law  from  the 
heats  of  formation  or  combustion,  available  for  many  compounds 
in  tabular  form.  Heat  capacities  of  most  simple  molecules  have 
been  calculated  statistically  to  very  high  temperatures  by  many 
authors  and  in  Table  7  are  given  energy  Integrals  derived  from 
such  data  in  analytical  form  accurate  to  about  Ifo  over  the 
temperature  range  2000  to  5000°K. 

Since  calculations  of  this  section  are  concerned  with  an 
"ideal"  state,  its  definition  is  arbitrary  and  will  be  chosen 
so  as  to  provide  the  greatest  convenience  in  calculations  and  yet 
to  facilitate  the  subsequent  problem  of  calculating  true  detona¬ 
tion  rates.  For  this  purpose  the  following  assumptions  are  made: 

(1)  All  product  gases  are  ideal  and  the  volume  of  solid  or 
liquid  products  may  be  neglected. 

(2)  All  dissociations  into  fre^e  radicals  and  atoms  are  non¬ 
existent. 

(3)  Only  solid  carbon,  CO,  COg,  H2O,  H2^02,N2  (and  gaseous 
sulfur,  SOg,  metals  or  metallic  oxides,  insofar  as  the 
corresponding  elements  are  present  In  the  explosive)  are 
formed . 

(4)  Oxygen  reacts  quantitatively  with  carbon  to  form  CO,  the 
excess  forms  quantitatively  H2O  and  what  is  left  over  reacts 
to  form  CO2. 

The-  question  of  whether  a  real  pressure  and  temperature  range 
exists  in  which  all  these  assumptions  hold  is  entirely  immaterial 
for  the  following,  because  and  other  "ideal"  quantities 
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resulting  from  these  calculat5.ons  are  merely  convenient  steps 
to  reach  the  real  detonation  rates. 

In  Appendix  II  of  this  report  examples  are  given  of  the 
detailed  procedure  followed  in  the  calculation  of  ideal  quantities 
here  defined, 

2 o _ The  calculation  of  real  detonation  rates . 

A  rigorous  calculation  of  true  detonation  rates  cannot  he 
carried  out  at  present,  a  sufficient  reason  for  this  being  that 
the  equation  of  state  of  gaseous  mixtures  at  densities  above  1  g/cc 
and  temperatures  of  over  5000*^11  is  unknown . 

In  fact  several  approximations  and  simplifications  and  one 
frank  assumption  (the  form  of  the  equation  of  state)  must  be  made 
and  they  will  now  be  considered  one  by  one. 

Dissociation  into  atoms  and  free  radicals  is  fortunately  not 
of  Importance  because  of  the  high  pressures  in  detonation  waves 
of  solid  explosives  at  ordinary  densities  of  loading.  Consider 
for  example  the  dissociation  H2  =  2  H.  At  5000OK  (rather  high 
for  most  explosives)  the  dissociation  constant  has  been  calculated 
statistically  to  be  44.7  Atm.  The  concentration  of  free  hydrogen 
does  not  exceed  10^  by  volume  with  most  explosives  and  if  the 
total  pressure  is  10^  Atm.  it  is  readily  found  that  2fo  of  hydrogen 
is  dissociated  into  atoms.  This  means  an  absorption  of  heat 
roughly  equal  to  5  Kcal  per  Kg  of  explosive,  whose  total  heat  of 
explosion  is  normally  of  the  order  of  1000  Kcal/Kg.  Thus  T2 
is  lowered  by  the  inclusion  of  this  dissociation  to  the  extent 
of  0.3^  but  at  the  same  time  the  mole  number  is  increased  and 
hence  the  pressure  is  greater.  This  effect  is  of  the  order  of 
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O.lOf;  and  hence  the  total  effect  on  the  detonation  rate  Is  (be¬ 
cause  of  the  square  root  relation)  (1/2)  (0.3-0.10)^&  =  0.10/. 
Similar  calculations  shov  that  the  water  dissociations: 

2H2O  =  Hg  +  2  0E  or  HgO  +  ,(l/2)02  =  20H  and 

2H2O  =  Og  +  2H2  are  of  small  importance.  Dissociations  of 

oxygen  and  nitrogen  are  still  more  insignificant. 

Various  polyatomic  molecules  have  been  reported  in  the  gas 
samples  withdrawn  from  bombs  after  detonation.  Thus  Schmidt 
(Part  I)  reports  CHi|,C2H2^  HCN,  C2N2,  in  variable  but  small 

amounts  from  almost  all  explosives.  To  allow  for  equilibrium 
formation  of  all  these  molecules  seems  mathematically  an  almost 
impossible  task,  particularly  when  it  is  considered  that,  the 
thermodynamic  functions  of  most  of  these  molecules  are  not  too 
accurately  known  at  the  high  temperatures  in  question.  Approximate 

I 

calculations  made  on  some  of  those  which  are  found  in  greatest 
amounts  indicate  that  the  observed  concentrations  are  in  excess 
of  equilibrium  under  the  conditions  of  the  detonation  wave  (thus 
with  TNT  Schmidt  reports  roughly  equimolecular  quantities  of 
ammonia  and  hydrogen  whereas  the  equilibrium  ratio  is  Issr  than 
IslO),  In  all  probability  these  molecules  are  formed  during  the 
cooling  process  after  detonation.  But  even  when  it  is  supposed 
that  the  analytical  data  represent  the  true  conditions  in  the 
detonation  wave,  the  formation  of  all  these  complex  molecules 
does  not  alter  greatly  the  results  of  calculation. 

Taking  as  an  example  the  particularly  unfavorable  case  of 
TNT  for  which  the  oxygen  deficiency  is  large  and  the  complex 
products  particularly  abundant  according  to  Schmidt,  one  finds 
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that  In  the  formation  of  all  complex  products  reported  by  Schmidt 
a  total  of  8.7  Ecal  of  heat  is  absorbed  per  Kg.  of  TNT,  vhili.e  the 
total  heat  evolution  is  375  Kcal,  At  the  same  time  the  mole 
number  of  gaseous  products  is  reduced  from  32.6  to  30*7  and  there¬ 
fore  the  total  effect  on  the  detonation  rate  is  approximately 
[ (30«7/32. 6) (875/866,3) ] =  0.965.  The  neglect  of  the  complex 
products  thus  overestimates  the  rate  by  3*3^,»  "bnt  it  is  believed 
that  in  general  the  corresponding  error  is  smaller. 

The  formation  of  free  hydrogen  in  the  detonation  vave  in  the 
absence  of  solid  carbon  can  be  neglected  because  the  equilibrium 
constant  of  the  reaction  HgO  +  CO  =  CO^  +  ranges  from  0,2 
at  3000°  to  0.07  at  5000°K.  Thus  the  greater  part  of  the  hydrogen 
is  present  as  vater  and  the  formation  of  a  few  moles  of  free 
hydrogen  per  Kg  of  explosive  is  of  little  effect  on  the  detonation 
rate.  This  follows  because  there  is  no  mole  number  change  in  the 
reaction  and  the  heat  evolution  is  only  4.5  Kcal,  per  mole  at 
3000°  and  is  less  at  5000*^.  The  neglect  of  hydrogen  in  the 
absence  of  solid,  carbon  causes  therefore  an  underestimation  of 
the  rate  by  a  few  percent  at  most. 

In  the  presence  of  free  carbon,  on  the  other  hand,  the  forma¬ 
tion  of  free  hydrogen  will  proceed  almost  quantitatively.  The 
equilibrium  constant  of  the  reaction  C  +  E2O  =  CO  +  H2  is  about 
105  Atm.  at  3000°K  and  may  be  estimated  as  5  x  10^  Atm,  at  5000°. 

Since  the  pressure  (or  fugacity,  to  be  more  correct)  of  carbon 

Id 

monoxide  in  the  detonation  wave  seldom  rises  to  10  Atm, ,  the 
ratio  /Pj^  0  greater  than  unity  so  long  as  carbon  is  present. 
The  assumption  that  the  reaction  goes  completely  to  the  right  means 
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an  overestimation  of  the  detonation  rate  because  the  heat  absorp¬ 
tion  is  overcompensated  by  an  increased  mole  number  of  gaseous 
products.  The  folloving  examples  in  Table  1  demonstrate  the 
niagnitude  of  the  effect  on  D  of  the  extreme  assumptions  concern-^ 
Ing  this  equilibrium.  It  is  seen  that  ideal  rates  calculated  with 
no  carbon  and  those  calculated  with  no  hydrogen  differ  only  by  3^ 
or  less.  In  actuality  the  first  case  is  approached  closely  almost 
always  and  hence  the  error  resulting  from  the  extreme  assumption 
must  be  small. 


Table  I 

Effect  of  equilibria  on  ideal  detonation  rates. 


Compound 

Decomposition  equation 

Ti°K 

Di  : 

Picric 

Acid 

6  CO  +  HgO  +  iH2  +  li 

5440 

2140 

5i  CO  +  li  E2O  +  |C  -r  1-1  Ng 

3615 

2120 

TWT 

6  CO  +  ^  Hg  +  C  +  l|  Mg 

2880 

2070 

CO  +  2  ^HgO  +  C  +  li 

Mg 

5650 

2005 

Tetryl 

7  CO  +  H2O  +  li  K2  +  2i  N2 

3950 

2560 

5i  CO  +  2|  HgO  +  li  C  +  2i 

OJ 

4245 

2290 

Summing  up  the  results  of  this  discussion  it  is  found  that 
the  effects  of  chemical  equilibria  in  the  detonation  wave  are 
divers  and  altogether  may  effect  a  decrease  of  the  rate  with 
increasing  pressure  (or  density  of  loading)  amounting  to  perhaps 
10  or  at  the  very  most  15^»  This  conclusion  is  interesting  from 
the  theoretical  point  of  view  since  it  indicates  that  the  study  of 
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detonation  rates  does  not  provide  strong  evidence  as  to  vhether 
chemical  equilibria  are  established  in  the  vave  or  not.  For  the 
detonation  to  comply  with  the  present  theory  it  is  only  necessary 
to  postulate  that  the  reaction  is  "complete"  in  the  sense  that 
the  oxygen  combines  vlth  all  the  carbon  to  form  CO,  the  excess 
oxygen  combines  with  hydrogen  to  form  water  and  if  any  is  still 
left  over,  it  combines  xrlth  CO  to  form  COg* 

Experimental  evidence  shows  that,  with  increasing  density  of 
the  explosive,  the  detonation  rates  rise  sharply  and  reach  values  al 
most  four  times  as  high  as  the  "ideal"  values  calculated  by  the 
described  procedure.  Other  simplifications  having  been  novr 
eliminated,  it  appears  probable  that  the  rise  is  due  to  deviations 
of  gases  from  ideal  behaviour,  unless,  of  course,  the  whole  theory 
is  fundamentally  Incorrect.  In  Part  I  the  calculations  of  Schmidt, 
who  has  made  use  of  the  Abel  equation  of  state,  have  been  reviewed 
to  show  that  qualitatively  at  least  deviations  from  ideal  behaviour 
may  be  made  to  account  for  the  observed  rise. 

However  the  Abel  equation,-  as  well  as  its  parent  the  van  der 
Vaals  equation  of  state,  is  specifically  derived  for  gases  at  low 
densities  and  must  fail  under  conditions  encountered  in  the  detona¬ 
tion  waves.  Bridgman  ha.s  studied  the  behaviour  of  nitrogen  at 
room  temperature  and  at  densities  as  high  as  1.2  g/cc  and  has  found 
that  the  semi-emplrlcal  equation  of  Becker  describes  his  experimenta 
observations  very  well.  This  equation  has  four  arbitrary  constants 
and  is  written  by  Bridgman  in  the  form: 

P=  (15) 
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At  the  high  densities  of  Bridgman’s  experiments,  however,  the  last 
two  terras  are  of  minor  significance.  Omitting  them  and  changing 
somewhat  the  value  of  the  constant  k  -  the  only  remaining 
arbitrary  constant  -  one  can  obtain  almost  as  good  an  agreement 
with  Bridgman's  data  as  he  was  able  to  obtain  with  four  constants. 
(The  second  correction  term,  describing  an  attraction  between 
molecules,  must  become  quite  unimportant  at  the  high  temperatures 
and  densities  of  the  detonation  waves  and  the  last  term  was  used 
by  Becker  to  obtain  agreement  at  the  critical  point.)  It  was 
decided,  therefore,  to  use  the  Becker  equation  in  abridged  form 
with  only  one  adjustable  constant.  It  is  known  from  several 
experimental  and  theoretical  sources,  however,  that  the  repulsive 
energy  between  molecules  can  be  represented  as  proportional  to 
the  inverse  ninth  to  twelfth  power  of  the  distance  of  approach. 
Now,  with  increasing  temperature  the  kinetic  energy  of  molecules 
Increases  linearly  and  hence  they  must  approach  closer  on  colli¬ 
sions.'  The  constant  k  of  the  Becker  equation  which,  as  one  can 
readily  convince  oneself,  becomes  equivalent  to  the  van  der  Vaals 
(or  Abel)  covolume  in  the  limit  of  low  densities  and  is  therefore 
proportional  to  the  cube  of  the  molecular  diameter,  must  be 

temperature  dependent.  Depending  on  which  of  the  molecular 

-1/3  -l/it 

repulsion  laws  is  chosen,  one  must  take  either  a  T  '  or  T  ^ 
dependence  of  k.  The  former  was  actually  chosen  for  the  calcula¬ 
tions,  because  it  fortuitously  shortens  the  algebraic  labor,  and 
the  final  form  of  the  equation  of  state  is:  (where  k  =  K/T  ) 
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In  the  form  In  which  the  equation  is  written  above,  the 
~l/5 

terra  KT  'is  the  covolurae  of  all  the  molecules  formed  from  M 
grams  of  the  reactant.  As  the  following  will  show,  it  may  be 
represented  as  a  summation  over  the  Individual  covolumes  of  the 
different  products  formed. 

The  substitution  of  equation  (l6)  into  the  basic  equations 
derived  In  Section  1  is  laborious  but  straightforward.  A  few 
simplifications  In  minor  terms  had  to  be  made, but  all  of  them  have 
been  checked  and  found  to  alter  the  result  insignificantly.  The 
calculation  is  discussed  in  detail  in  the  Appendix  and  here  only 
the  final  results  for  the  ratio  D/D^  and  for  the  equivalent  of 
equation  (8)  will  be  written 


In  these  equations  the  variable  x  has  been  used  to  abbreviate 
-1/3 

KT  /vg,  while  »»  -  '?k/2X Since  equation  (l8)  re¬ 

lates  the  quantity  k^,^'/jto  x  ,  the  ratio  D/D^^  may  be  calculated 
and  plotted  as  the  function  of  the  former.  In  making  this  calcula- 
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tion  it  is  found  that  the  numerator  of  the  first  factor  3.n  eque- 
tion(33/  contributes  by  far  the  greatest  variation  in  D,/tjL  vith 
density.  The  product  of  all  the  other  factors  makes  at  most 
difference  in  the  result.  In  equation  (3^)  the  effect  of  reason¬ 
able  changes  in  the  values  of  ^  and  (^/^t  )  is  extremely 
small.  This  means  that  the  actual  values  of  the  heat  capacity 
of  the  product  mixture  are  of  very  small  import  for  the  shape  of 
the  D/D^  curve.  In  Table  2  three  calculations  have  been  made, 
the  first  for  the  second  for  =  9  and  the  third 

for  =  6,5  cal/mole  degree.  All  knovrn  explosives  fall  into 

this  range  of  Cyj_  values  and  the  variation  between  the  three  cases 
is  so  small  that  if  no  great  accuracy  were  desired  one  could  use  a 
universal  function  for  all  explosives.  In  Table  3  ve  give  the 
values  of  kij  and  T/T^.  They  depend  only  very 

slightly  upon  the  value  of  chosen  (the  calculations  have  been 
made  for  =  8)  and  therefore  the  standard  values  of  {X/^\)  were 
used  for  the  second  factor  in  equation  (35 )•  A  further  simplifica¬ 
tion  in  the  calculations  was  the  neglect  of  the  dependence  of  k, 
and  upon  T  as  it  varies  with  increasing  density  of  the 
explosive.  The  corresponding  changes  are  very  small  and  so  long 
as  a  consistent  procedure  is  used  for  all  explosives  the  damage 
to  the  results  is  entirely  insignificant. 

In  Table  4  have  been  collected  practically  all  published 
determinations  of  the  detonation  velocities  of  militarily  interes¬ 
ting  pure  and  mixed  explosives.  Only  measurements  on  ’’safety'' 
explosives  have  been  excluded  as  these  are  usually  made  up  of 
coarse  particles  of  the  combustible  and  the  oxygen  carrier,  such 
inhomogeneity  being  deliberately  resorted  to  in  order  to  slow  down 
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the  detonation  velocity  belov  that  demanded  by  the  theory.  For 
each  explosive  the  Ideal  temperature  and  detonation  velocity  vere 
calculated  by  the  method  outlined  In  Appendix  II  and  the  ratios 
D/D^  formed  using  experimental  data.  Then  that  one  of  the  three 
functions  given  In  Table  2  was  used,  whose  Cvi  value  was  nearest 
to  the  calculated  In  the  detonation  wave;  from  It  the  values 
of  k  /M  were  read  off  for  all  experimental  data.  The  correspond¬ 
ing  k  are  then  easily  calculated  and  these  are  given  In  the  3  rd 
column  of  the  Table.  The  averages  of  the  k  for  each  explosive 
are  also  given;  the  deviations  of  individual  values  from  such 
averages  are  rather  large.  It  must  be  considered,  however,  that 
experimental  values  of  both  the  detonation  velocity  and  the  bulk 
density  enter  Into  the  calculation  and  in  view  of  the  nature  of 
such  experiments  accidental  errors  of  quite  a  few  percent  are  to 
be  expected.  This  Is  supported  by  the  observation  that  the  data 
of  some  Investigators  are  uniformly  much  better  than  those  of  the 
others  and  that  for  more  thoroughly  studied  explosives  better 
agreement  is  obtained. 

But  over  the  accidental  scattering  there  is  superimposed  a 
very  definite  trend  of  k  which  decreases  with  increasing  density 
of  loading.  This  suggests  that  the  equation  of  state  taken  as  the 
basis  of  these  calculations  Is  not  exactly  valid  over  the  pressure 
range  covered  in  the  experiments.  A  better  agreement  might  have 
been  obtained  by  the  use  of  an  equation  of  state  requiring  a 
slightly  greater  dependence  of  covolume  upon  pressure.  However, 
the  available  experimental  detonation  rate  data  are  not  accurate 
enough  to  justify  a  further  elaboration  of  the  theory  at  the 
present  time.  V/lth  the  average  values  of  k  given  for  each 
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explosive  one  can  obtain  a  representation  of  the  detonation  velocity 
curve  vhich  is  vlthln  5^  of  the  mean  of  experimental  data  over  the 
density  range  of  0.7  to  1.7  g/cc  for  organic  explosives. 

This  result  must  undoubtedly  be  regarded  as  a  considerable 
triumph  of  the  theory  and  from  the  theoretical  point  of  view  it  is 
valuable,  being  the  first  Information  on  the  behaviour  of  gases 
under  conditions  which  might  best  be  termed  as  those  half-way  be¬ 
tween  conditions  realisable  in  the  laboratory  and  those  existing 
in  the  stars. 

Prom  the  practical  point  of  view,  however,  the  theory  is  of 
little  Interest  unless  it  provides  the  means  for  calculating  de¬ 
tonation  velocities  for  new  compounds.  Fortunately  the  theory  is 
capable  of  accomplishing  this  task  as  the  following  will  con¬ 
clusively  show. 

In  Table  5  have  been  collected  the  average  values  of  k  from 
Table  4,  the  second  column  giving  these  values  corrected  to 
with  the  aid  of  the  law.  The  third  column  gives  the  cal¬ 

culated  values  of  k^QQ  obtained  from  a  set  of  covoluraes  for  the 
individual  products  given  in  Table  6.  These  values  have  been  ob¬ 
tained  by  a  process  of  gradual  fitting  and  could  be  slightly 
altered  without  Impairing  the  excellent  degree  of  agreement  with 
experimental  data  which  is  shown  in  Table  5.  Since  the  covolumes 
of  CO,  N2  and  O2  vere  chosen  as  identical,  there  remained  four 
arbitrary  parameters  to  adjust  for  the  best  possible  fit.  The 
number  of  compounds  for  which  excellent  agreement  with  experimental 
data  is  obtained  is  greater  than  ten  and  the  degree  of  agreement 
is  as  good  as  the  quality  of  experimental  data  allows.  Therefore 
the  agreement  is  hardly  an  acclrJental  one  and  full  confidence  may 
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be  placed  in  calculations  of  k's  of  nev  organic  compounds  vlth 
this  set  of  individual  covoluraes  unless  solid  carbon  is  formed  be¬ 
cause  of  oxygen  deficiency.  Some  cases  of  explosive  mixtures  have 
been  included  in  the  Tables  and  for  these  an  equally  good  agreement 
of  theory  and  experiment  is  obtained.  Thus  one  may  s&telj  apply  the 
present  theory  to  mixtures  of  explosives  vhich  might  be  developed 
in  the  search  for  powerful  and  yet  safe  compositions. 

The  individual  covolumes  given  in  the  Table  show  a  surprisingly 
good  agreement  with  the  van  der  ¥aals  constants  of  these  gases 
which  have  been  obtained  from  critical  data  and  by  other  methods 
(the  critical  data  constants  were  all  recomputed  to  300°K  using 
the  same  law  as  for  the  detonation  data).  The  detonation 

data  seem  to  be  somewhat  higher  on  the  average,  but  this  in  no  way 
detracts  from  the  success  of  the  agreement,  because  extrapolation 
over  a  15  fold  range  of  temperature  had  to  be  made  in  arriving  at 
the  300*^  data  and  a  choice  of  only  a  sllghtlj’’  different  exponent 
of  T  (for  Instance  T  )  might  have  eliminated  the  difference. 

For  steam,  however,  the  selected  value  is  relatively  too  small;  on 
the  other  hand  it  is  well  known  that  the  law  of  corresponding 
states  holds  less  well  for  steam  than  for  other  molecules  here 
considered.  The  change  from  critical  data  for  steam  is  actually 
in  the  right  direction.  For  carbon  dioxide  on  the  other  hand  an 
excessively  large  covolume  had  to  be  taken  while  for  h^^'drogen  it 
is  too  small  and  we  ca.n  think  of  no  good  reason  why  this  should  be 
so.  For  solid  carbon  l4  cc/raole  was  taken  because  at  4000°K  this 
becomes  about  6,  thus  giving  the  experimental  density  of  graphite 
(2.2  at  room  temperature)  under  the  conditions  in  the  detonation 
-’-ave.  In  the  case  of  solid  carbon  (and  liquid  lead,  see  below). 
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the  use  of  the  expression  covolume  is  really  improper.  To  the 
extent  that  the  gases  are  insoluble  in  the  condensed  phase  it  is 
the  real  volume  of  the  latter  per  gram  mole  vhich  enters  into  the 
calculations  and  which  is  given  in  Table  6.  The  assumption  of 
the  temperature  dependence  is  obviously  erroneous  and  no 

physical  significance  should  be  attached  to  values  given  in  Table  6 
for  300°K,  These  are  there  merely  for  the  purpose  of  uniformity. 

The  data  for  TNT  are  definitely  out  of  line  with  the  rest  of 
the  experimental  figures  whether  one  chooses  hydrogen  or  water 
vapor  as  the  Ingredient  of  the  detonation  wave.  It  is  to  be  noted 
however  that  these  extreme  assumptions  brachet  the  k  calculated 
from  individual  covolumes.  Fovr,  TNT  is  a  rather  cool  explosive 
with  a  temperature  of  about  3000°  in  the  detonation  wave  and  at 
this  temperature  the  equilibrium  constant  of  the  reaction 
C  +  H2O  =  CO  +  H2  is  such  (c.f.  an  earlier  discussion)  that  at 
the  high  pressures  of  the  detonation  wave  only  partial  conversion 
according  to  the  above  equation  is  to  be  ejcpected.  The  extent  of 
this  conversion  will  vary  of  course  with  the  density  of  the 
explosive  and  it  is  very  difficult  to  estimate  because  of  the 
non-ideality  of  the  gas  mixture.  The  partial  conversion  has  little 
effect  on  the  detonation  velocity  as  was  seen  in  Table  1,  but  the 
experimental  value  for  k  should  lie  between  the  values  calculated 
from  the  individual  covolume  constants  for  the  extreme  cases  con¬ 
sidered,  This  is  as  found  experimentally  and  although  the  treat¬ 
ment  is  not  quantitative  and  hence  cannot  be  regarded  as  particularly 
successful,  it  shows  nevertheless  that  TNT  is  not  in  disagreement 
with  the  theory.  It  is,  however,  a  good  illustration  of  the  remark 
nade  in  the  beginning  of  this  report  that  when  one  of  the 
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ingredients  of  the  detonation  wave  is  divided  between  two  phases, 
the  theory  ceases  to  be  very  effective.  With  all  other  explosives 
tried,  except  possibly  TKB,  enough  oxygen  is  available  to  oxidize 
all  carbon  to  carbon  monoxide  and  at  least  some  hydrogen  to  water. 
The  ratio  is  then  such  that  no  reduction  to  solid 

carbon  can  take  place  even  at  J>000°  or  slightly  below.  Solid 
carbon  found  analytically  among  the  detonation  products  of  these 
explosives  must  be  formed  in  the  cooling  process. 

The  treatment  of  lead  azide  leaves  also  something  to  be 
desired.  One  possible  assumption  is  that  all  lead  is  vaporised, 
in  which  case  one  obtains  Tj_=5300°K  and  Dj_=1250  m/sec.  The  boil¬ 
ing  point  of  lead  is  about  1900°K  and  thus  at  3300°  it  should  have 
a  vapor  pressure  of  only  about  100  Atm.,  quite  insignificant 
compared  with  the  total  pressure  of  the  order  of  lO-^Atm.  This 
assumption  is  definitely  erroneous  therefore  and  our  attempts  to 
Improve  it  by  supposing  that  lead  vapor  is  diatomic  did  not  change 


this  conclusion.  One  must  then  assume  liquid  lead,  which  gives  a 
much  higher  temperature  in  the  detonation  vrave  (because  of 
greater  heat  evolution)  T^=  4600°K,  but  substantially  the  same 
detonation  rat,e,  Dj_=1275  because  of  reduced  mole  number  of  gaseous 
products.  The  calculation  was  carried  out  assuming  for 

liquid  lead  of  6  cal  per  gram  atom  and  a  heat  of  fusion  of  2  Kcal/g. 
atom  but  a  reasonable  change  in  these  assumptions  will  not  alter 
the  by  more  than  a  very  few  percent.  The  vapor  pressure  of  lead 
under  those  conditions  is  undoubtedly  much  higher  than  the  100  Atm. 
calculated  for  3300°,  but  in  all  likelihood  it  is  still  insignificant 


compared  with  the  total  pressure  in  the  wave.  Partial  evaporation. 


In  any  case,  does  not  seriously  alter  the  value  of  Dj_,  since  it 
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results  in  a  decreasing  temperature  but  an  increasing  mole  number 
of  gaseous  products  and  these  happen  to  compensate  each  other 
almost  quantitatively  in  the  expression  for  Dj_.  Prom  the  ex¬ 
perimental  k  one  finds  a  volume  for  liquid  lead  of  115  cc/gr.atom 
at  or  about  45  cc  at  the  temperature  of  the  detonation  wave. 

This  means  a  density  of  lead  of  about  4.5>  a  value  which  might 
appear  quite  unreasonable  but  really  is  not.  Thus  for  mercury  (c.f. 
Handbook  of  Chemistry,  Chem,  Rubber  Co.)  a  density  between  4  and  5 
is  reported  for  the  critical  temperature  and  4600°K  should  be  not 
far  from  the  critical  temperature  of  lead.  The  result  of  these 
calculations  is  therefore  an  entirely  reasonable  one,  but  of  course 
it  is  very  doubtful  whether  the  theory  in  its  present  form  can 
justifiably  be  applied  to  liquids  near  their  critical  temperatures 
and  caution  must  be  advised  for  all  other  explosives  that  may 
yield  large  amounts  of  liquid  products  in  the  detonation  wave, 

4,  The  search  for  more  powerful  explosives. 

It  seemed  to  us  to  be  of  some  interest  to  establish  rules 
that  should  govern  the  synthesis  of  new  and  more  powerful  explosives 
but  unfortunately  the  theory  is  so  complex  that  it  is  Impossible  to 
state  such  rules  in  a  simple  explicit  form. 

To  obtain  the  highest  possible  one  must  have  a  compound 

with  the  greatest  possible  heat  of  formation  (as  endothermic  as 
possible)  and  one  that  gives  the  greatest  possible  number  of  moles 
of  products  with  the  largest  possible  heat  evolution.  From  the 
considerable  amount  of  calculation  that  has  been  done  by  us  it 
appears  that  the  oxygen  balance  need  not  be  so  high  as  to  give 
mainly  CO2  :  the  greater  heat  evolution  in  forming  this  compound 
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is  compensated  by  the  smaller  number  of  moles  formed  per  unit  weight 
of  the  explosive  and  by  the  higher  heat  capacity  of  carbon  dioxide. 

One  of  the  highest  thus  far  found  is  that  of  Cyclonite  from 

which  no  carbon  dioxide  is  formed  at  all.  On  the  other  hand,  when 
the  oxygen  balance  sinks  below  the  amount  needed  to  give  complete 
oxidation  to  CO  and  H2O,  then  the  temperature  in  the  detonation  be¬ 
comes  so  low  that  small  values  of  are  the  inevitable  result. 

When  it  comes  to  the  calculation  of  real  detonation  rates,  the 
prediction  becomes  much  more  difficult  because  additional  factors 
come  into  play.  Steam  and  hydrogen  have  very  small  covolumes  and 
therefore  give  small  D/D^  ratios;  COg  on  the  other  hand  has  a  very 
large  covolume  and  is  therefore  favorable.  High  temperatures  in 
the  detonation  wave,  which  give  large  also  reduce  the  magnitude 
of  the  covolumes  and  so  result  in  smaller  D/Djj^  ratios.  To  give 
hard  and  fast  rules  to  follow  in  the  search  for  the  ideal  explosive 
under  these  circumstances  would  be  presumptuous.  It  is  suggested 
instead  that  a  calculation  following  instructions  given  in  Appendix  II 
be  made  after  the  conception  of  a  new  explosive  on  paper  and  before 
its  synthesis  in  the  laboratory. 


Table  2 
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Eela1>ion  between  D/d^^  and  S  f  l/M 


k  (pi/M 

Case  I 

(Cy  =  11..5) 

Use  when  ^£•1.22 

Case  II 

(Ov  =  ?)  — 

Use  when  1.28  ^1.22 

Case  III 
(Gy  =  6.5) 
Use  when  -  1 

0 

1 

1 

1 

0.05 

1.10 

1.09 

1.08 

0.10 

1.19 

1.18 

1.17 

0.15 

1.29 

1.28 

1.26 

0.20 

1.39 

1.37 

1.35 

0.25 

1.49 

1.47 

1.45 

0.30 

1.60 

1.58 

1.55 

0.35 

1.71 

1.69 

1.66 

o.4o 

1,84 

1.81 

1.77 

0.45 

1.96 

1.92 

1.88 

0.50 

2.09 

2.05 

2.00 

0.55 

2,22 

2.17 

2.12 

0.60 

2.35 

2.30 

2.24 

0.65 

2.49 

2.43 

2.36 

0.70 

2.64 

2.57 

2.50 

0.75 

2.79 

2.71 

2.64 

0.80 

2.54 

2.86 

2.78 

0.65 

3.10 

3.01 

2.92 

0.90 

3.27 

3.17 

3.07 

0.95 

3.43 

3.34 

3.22 

1.00 

3.61 

5.51 

5.39 

1.05 

3.79 

5.68 

5.55 

1.10 

3.98 

5.86 

5.71 

1.15 

4.18 

4.05 

3.88 

1.20 

4,39 

4.25 

4.06 
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Use  linear  interpolation  between  the  values  given. 


Table  3 
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!3C 

Variation  of  T/^j  ,  of  Tg/^  and  of  a  function  (x  +  x  e  ) 
with  Loading  Density  as  Calculated  on  the  Basis  of  Gas  Im¬ 
perfections. 


^,k/M 

^2/^21 

2  ; 

X  +  X  e 

0.000 

1.000 

1.000 

0.0 

0.100 

0.971 

0.980 

0.4 

0.200 

0.948 

0.962 

0.9 

0.500 

0.930 

0.944 

1.6 

0.400 

0.912 

0.925 

2.4 

0.500 

0.899 

0.907 

3.8 

0.600 

0.886 

0.889 

5.6 

0.700 

0.874 

0.871 

7.6 

0.800 

0.864 

0.855 

10.0 

0.900 

0.854 

0.855 

12.8 

1.000 

0.845 

0.816 

16.5 

1.100 

0.836 

0.798 

20.6 

1.200 

0.829 

0.780 

26.0 

Table  4 
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Covolumss  from  experimental  data. 


Density  Detonation  Experimental  Deviation  from  Mean 

rate  ^  mean,^ 

1  o 

EETE.  M  «  516^  Hfo  =  -123,000  cal/mole*  T^  =  ^930  K;  Dj_  =  2390  m/sec. 
•  Case  I 


Data  of  Friedrich 


0.50 

3960 

204- 

10 

' 

0.65 

1+1+00 

195 

5 

0.80 

1+900 

I9I+ 

5 

1.00 

5500 

184 

-1 

1.20 

6300 

l£5 

-1 

1.40 

7100 

1® 

-1 

1.60 

7900 

181 

-2 

Data  of  Both 

0.1+5 

3510 

172 

7 

185  cc  at  493 0°K 

0.75 

1+820 

196 

6 

0.91 

5300 

195 

4 

1.01+ 

5730 

188 

2 

1.1+5 

7375 

184 

-1 

1.72 

8050 

171 

-7 

Data  of  Tonegutti 

1.20 

6070 

176 

-5 

1.56 

7630 

178 

.4 

1.60 

8200 

Data  of  others 

IH7 

1 

1.71 

85I+O 

178 

-4 

Cyolonite  (Hexogen)j  M  =  222;  HfQ  =  21,300; 

T* 

=  4650°K;  D^  =  2510 

Case  II 

Data  of  Tonegutti 

1.20 

6310 

126 

2 

1.56 

7890 

122 

-2 

124  cc  at  4650Ok 

Data  of  others 

1.67 

®80 

126 

2 

3 -nitroglycerine;  M  =  227;  Hf^  =  -82,700;  T 

i  = 

5230°;  Di  =  256  0 

Cass  I 

Data  of  Friedrich 

1.59 

8020 

- 133— 

2 

1.75 

solid,  8210 

126 

-4 

131  cc  at  5230®K 

Data  of  others 

1.59 

81+50 

Ul 

T 

1.59 

7I+50 

123 

-6 

1.  Heat  of  formation  from  the  elements. 


Oovolumes  from  experimental  data  (continued) 
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Tetryl;  M  =  287;  Hfo  «  9,500;  =  5950°;  =  2560 

Case  III 


Data 

of  Itiedrich 

1.00 

5480 

Ibl 

4 

1.28 

6500 

177 

2 

1.45 

7210 

176 

2 

1.54 

7380 

170 

-2 

1.61 

7480 

165 

-5 

Data 

of  Both 

0.61 

2570 

235  (omit) 

+56^ 

0.92 

4810 

159 

-8 

1.24 

6400 

178 

5 

1.55 

7000 

185 

7 

1.46 

7510 

1® 

6 

1.58 

7665 

172 

0 

Data 

of  others 

1.55 

7145 

165 

>4 

1.59 

7160 

159 

-8 

1.65 

7200 

154 

-11 

1.68 

7485 

158 

-8 

0.96 

5420 

185 

7 

Picric  Acid;  M  = 

229;  Tj  =  3440;  flfo  =  -55 

,500;  Di  = 

Cass  III 

Data 

of  Briedrich 

1.03 

5150 

146 

4 

1.23 

5820 

155 

-4 

1.59 

6450 

155 

-4 

1.63 

7210 

151 

‘7 

Data 

of  Both 

.18 

1.10 

4560 

115 

1.23 

5430 

146 

5 

1.34 

6110 

l4l 

0 

1.47 

6940 

149 

6 

1.47 

6060 

127 

-10 

1.56 

7500 

152 

8 

Data 

of  Kast 

1.34 

6160 

142 

1 

1.46 

6700 

144 

2 

1.55 

7000 

143 

1 

1.60 

7100 

159 

-2 

Data 

of  Dautriche 

0.90 

4900 

158 

12 

1.20 

5800 

147 

4 

1.50 

7100 

149 

6 

Data 

of  Forg 

0.85 

4629 

155 

8 

1.20 

5215 

129 

-8 

1.4o 

6232 

159 

-1 

1.60 

7021 

158 

-2 

Data 

of  others 

0.98 

4970 

147 

4 

1.65 

7600 

147 

4 

175  cc  at  3950°K 


2140 


l4l  cc  at  3250°K 


Covolumea  from  experimental  data  (continued) 
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TOT;  M  =  227;  Hfo  =  -15,000;  =  2950°K;  Dj  =  2110; 


(or  Tj 

=  5650°  ; 

D^  =  2005  if  water  is 

quantitatively  formed) 

Case  III 

Data 

of  Eriedrich 

1.00 

4700 

126 

-6 

1.29 

5900 

142 

6 

1.46 

6500 

140 

5 

1.57 

6749 

156 

1 

1.59 

6900 

158 

5 

Data 

of  Both 

154  cc  at  2950^^ 

1.21 

4720 

112 

-16 

1.27 

4650 

105 

-22  or 

1.40 

5900 

150 

-5  l4l  cc  at  5650°K 

if  H^O  is  formed 

1.52 

7400 

154 

15 

Data 

of  East 

1.54 

5940 

137 

2 

1,45 

6400 

153 

5 

1.50 

6590 

158 

5 

1.60 

6680 

132 

-2 

Data 

of  Dautrlche 

0.90 

4100 

121 

-10 

1.20 

5600 

145 

7 

1.40 

6600 

149 

11 

Data 

of  Forg 

0.84 

4400 

145 

8 

1.22 

5596 

154 

0 

1.45 

6294 

156 

2 

1.59 

6772 

155 

1 

Data 

of  others 

1.62 

7000 

137 

2 

1.66 

6565 

121 

-10 

Nltropentanone  (0^2.2^15%^? 

r  584;  Hfo 

=  -155,000;  =  5350;  D^  =  2250 

Case  III 

Data 

of  Wohler 

1.15 

4650 

180  (omit) 

-26 

1.50 

6020 

22S 

-6 

1.44 

7170 

250 

1  244  cc  at  5550°K 

1.59 

7940 

254 

4 

Nitropentanol  (0^2.5  *^15^5  ^  = 

‘‘31;  Ef„  : 

=  -140,000;  Tj  c  4090°;  Dj^  =  2400 

Case  II 

Data 

of  Wohler 

0.75 

5060 

510  (emit) 

26 

0.91 

5100 

258 

5 

1.01 

5800 

279 

1.11 

5940 

258 

.  0. 

5  246  cc  at  4090°K 

1.29 

6100 

229 

-7 

1,51 

7050 

25^ 

-5 

...57 

7560 

256 

-4 

;..62 

7540 

257 

-4 

Covo].ixmes  from  expsrimsntal  data  (continued) 
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Nitrohexanone  )|  M  =  598;  =  -163,000; 

Case  III 

Data  of  Wohler 

l.l8  6000  '  269  3 

1.25  5710  2J+2  -8 

1.1+2  7000  269  2 

1.51  77'+0  279  6 

Nitrcacenanol  (Cio%5°15^5^»  ^  ^  ="152,000,-  T, 


Case  III 

Data  of  Wohler 


0.81 

5470 

535  (omit) 

16 

1.00 

5820 

299 

4 

1.28 

6800 

285 

-1 

1.44 

7670 

289 

0 

1.52 

7900 

283 

-2 

TWB;  M  = 

=  213;  Hf^ 

,  =  -2,300;  Ti  =  5550;  Dj 

=  2260 

Case  III 

Data 

of  Dautriche 

0.80 

i+100 

112 

-10 

1.10 

5400 

128 

+5 

1.55 

6700 

158 

+  12 

Data 

of  others 

1.60 

7000 

122 

-2 

1.65 

7000 

119 

-4 

Mercury 

fulminate 

1;  M  =  284;  Hfo  =  65,000; 

Ti  =  5550' 

Cass  III 

Data 

of  Roth 

2.21 

5570 

74 

8 

2.70 

4000 

71 

4 

3.27 

4l60 

62 

-8 

3.71 

4810 

66 

-5 

Data 

of  others 

1.69 

3200 

80  (omit) 

17 

1.81 

3054 

69 

2 

1.82 

3122 

72 

4 

2.15 

5504 

67 

-1 

2.98 

4125 

67 

-0 

2.99 

5971 

63 

-7 

5.60 

4500 

62 

-8 

Lead 

azide;  M  =  292;  Hf^  = 

Data  of  Roth 

-105,000;  ST 

Cass  III 

4640  ;  = 

2.65 

3850 

101 

7 

2.75 

4420 

111 

19 

3.25 

4700 

99 

6 

3.50 

4850 

95 

2 

■T.  =  2880;  Di  =  2150 


265  cc  at  2880^ 

=  5820°;  D^  =  2I+OO 


289  cc  at  5820°K 


124  cc  at  5350^^ 


Di  =  1640 


67.5  cc  at  5530°K 
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Covolumes  from  experimental  data  (continued) 

Lead  azide  (continued) 

Data  of  Friedrich  (in  lead  tube  only;  glass  confinement  gives 
utterly  different  results;  seems  prolatle  that  shock  wave  through  glass 


rather 

than  the  detonation 

rate  in  the 

salt  was  ms 

jasured) . 

2.56 

4480 

118 

27 

3.51 

4745 

89 

.4 

4.05 

5275 

88 

-5 

^  cc  at  464o°K 

Data  of  Stadler 

1.45 

2120 

70 

-25 

deviations  from  the 

2.16 

3080 

90 

-3 

mean  are  here  much 

3.12 

4270 

93 

0 

greater  than  with  other 

3.19 

4540 

96 

5 

explosives  investigated 

Data  of  others 

3.44 

4500 

65 

-30 

Di-PETIT  (03122018%);  M  =  480;  Hf^  =  -210,000;  =  463  D^  =  2340 

Case  I 

Data  of  Friedrich 

1.52  7000  276  279  cc  at  463 0°K 

1.59  7410  282 


1.50 


Covolumes  frcm  experimental  data  (continued ) 

51 

-illo  EETN;  1Q$  K;  %  DNT;  (iKg);  Ef.  =  -704,000;  T.  =  4680°;  D.  =  2400 
485^  N%N05 

Case  I 

Data  of  Eaouai 

1.36  6600  542  oc  at  468o°K 


Table  5 
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Calculation  of 

ccvolumes 

per  mole  of 

explosive. 

Compound 

T^*^K 

k^,cc 

^00 

4500 

Bsmarks , 

Deviation, 

EETN 

Ij-930 

185 

1+72 

459 

-3 

Cyclonite 

1+650 

I2I+ 

310 

308 

-1 

Nitro¬ 

glycerine 

5250 

151 

559 

535 

-2 

Tetryl 

5950 

173 

1+08 

405 

-1 

Picric  Acid 

51+40 

ll+l 

318 

320 

+1 

Wet  P.A.* 

3200 

155 

340 

333 

-2 

TNT 

2950 

151+ 

287 

322 

If  no  HgO,  only  is 

formed 

TNT 

5650 

ll+l 

324 

If  only  H2O,  no  H2  is 

formed 

Nitropenta- 

none 

5550 

2I+J+ 

545 

543 

-1 

Nitropenta- 

nol 

1+090 

2l+6 

586 

604 

+5 

Nitrohexa- 

none 

2880 

265 

560 

573 

+2 

Nitrohexanol 

3820 

289 

672 

651 

-6 

TNB 

3550 

12I+ 

283 

302 

+7  Probably  some  car¬ 
bon  and  vater 

formed  as  with  TI-IT 

Hg(0CN)2 

5550 

67.5 

178 

- 

PbNg 

1+61+0 

95 

252 

Di-PSTN 

1+650 

279 

695 

705 

+1 

Glycol-di- 

nitrate* 

5500 

90 

254 

230 

-2 

1+970 

106 

270 

274 

+1 

Dry  guncotton*  5720 

613 

1420 

1456 

+2  for  1  kg. 

Wet  guncotton*  l880 

7I+2 

1370 

l4l0 

+3  ” 

PBTN-NH^NO^- 

mixture* 

NG-DNT 

1+680 

51+2 

1360 

l446 

+6  for  1  4g. 

*  Only  one  experimental  value  available  for  each  of  these. 
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Cavolumes  of  different  molecules. 


Substance 

0 

0 

0 

0 

van  der 
Waals 
at  c.p. 

corrected 
to  500° 

covolume  from 
viscosity 

Kr, 

C 

58.7 

58.6 

30.6 

.  39 

CO 

38.7 

39.9 

31.9 

39 

°2 

38.7 

31.8 

26.0 

51  —  enters  only 
into  nitro¬ 
glycerins 
data 

0 

0 

ro 

67.0 

42.6 

42.5 

44 

-SqO 

25.8 

30.5 

38.5 

21,2 

^2 

7.75 

26.6 

12.7 

25.7 

Hg  62  :  chosen  to  fit  Hg(0CK)2;  for  a  close  packing  of  spheres 

gives  a  density  of  about  9  instead  of  experimental  of 
about  l4  for  solid  mercury 


C  (solid)  l4  :  chosen  to  give  a  density  of  about  2  at  4000°K 

Pb  (liquid)  II6  ;  gives  a  density  of  about  4.5  at  4500°j  fits  PbKg  data 


TaUe  7 
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Empirioal  Equations  to  Bapresent  Values  of 
\  CvdT  ~  Cvi(l2~300)  in  the  Benge  2000°  5000°K. 


C_  =  fi  +  tTo  • 

Ti  2 


Substance  Equation  Accuracy 


C  (solid) 

-  -3 

=  5.9^  +  0.27  X  10  •^Tg 

1% 

00 

C  =  5.61  +  0.21  X  10“^T^ 

COg 

=  10.16  +  OM  X  lO'^Tg 

li 

^2 

_  -3 

=  4.91  +  0.30  X  10  Tp 
i 

HgO 

-  -3 

0^^  =  7.86  +  0.55  X  10  Tg 

li 

^2 

=  5.54  +  0.21  X  lO'^Tg 

Vfo 

Appendix  I. 


Calculation  of  the  ratio  P/P-t . 


A  division  of  equation  (7)  by  equation  (l4)  leads  to  an  expres¬ 
sion  for  the  ratio  D/D^: 

Va 


D  - _ 

'Si  "fKTi) 


(19) 


in  vhich  all  quantities  may  in  principle  he  calculated  from  equa¬ 
tions  (8,9#10,12  and  13) •  This  is  true,  of  course,  only  if  it  is 
assumed  that  chemical  equilibria  do  not  shift  the  composition  of 
the  explosion  products  with  density  of  the  explosive.  A  discussion 
given  in  Section  3  of  the  report  has  shown  that  the  effects  of  such 
equilibrium  shifts  upon  the  detonation  velocity  are  very  small  and 
it  has  been  decided  therefore  to  neglect  them  altogether  in  this 
calculation.  The  ratio  n/n^^  is  assumed  therefore  to  be  unity. 
Actually  it  is  somewhat  different  from  unity  (and  the  ratio  T/T^ 
is  not  what  we  calculate  it  below  due  to  the  same  cause)  but  these 
small  effects  are  automatically  taken  care  of  to  some  extent  by  the 
method  of  calculation,  which  blames  all  variations  in  on 

gas  imperfections.  By  comparing  equations  (5)  and  (16)  we  note 
that; 


Hence,  if 


-f  (V,  I )  = 


K 


=  X 


fV) 


(1  + 


(20) 


(21) 
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and: 


If  the  folloving  abbreviation  is  introduced: 

„  E _ r‘  r 

Cv  -  gyf'c'v)  ~  Ov  and  ~2.ll' ’ 


and  furthermore  a  superscript  bar  is  used  to  Indicate  the  averaging 

and 

of  C.^  over  the  temperature  range  Tp  *  'Tn  /components  of  the 
mixture,  then  (c.f.  equations  (9),  (lO)  and  (l^)  )  the  expression 
for  T2  may  be  written  thus: 

-r  H 

CC  r^CCf 

In  the  second  term  has  been  replaced  by  because  the 

second  term  is  much  smaller  tha.n  the  first  a.nd  the  total  varia¬ 
tion  '  from  C.^j_  to  a.raounts  to  only  some  15/  at  the  highest 

densities.  Prom  equation  (11 )  the  energy- volume  coefficient  is 
readily  evaluated: 
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It  remains  still  to  calculate 

;Cp= 


Now: 


(28) 


and  the  equation  of  state  gives:  . 

,  ,  I  ,  V,  XO*_  '/g 

(iii\  -  ^  ^  ^ 

VST/p  -  T  r+ 


\dl  /v/ 


tJ,  ^rSfe^-  -  “‘®'J  • 


(2^ 


(30) 


(31) 


These  expressions  may  he  substituted  into  the  equations  for  Up  and 
C„  and  thus  'P  obtained  as  a  function  of  C,,,  and  the  variable  x. 
This  value  of  may  in  turn  be  substituted  into  (2?)  so  that  the 

temperature  ratio  also  becomes  a  function  only  of  ideal  heat 
capacities  and  of  the  variable  x. 

Heturnlng  now  to  equation  (8)  and  using  the  derived  expressions 
for  f’(v)  and  f(v)  one  can  rewrite  it  thus: 


X 


I 


I  +- 


But  l/v^  =  ^  /M  and  if  both  sides  are  multiplied  by 

one  obtains  the  relation: 


(32) 


‘A. 


J<_  P- 


X 


(35) 
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vhere,  of  course,  the  complete  expression  for  should  he  sub¬ 

stituted.  Equation  (33)  is  the  same  as  equation  (l8)  of  the  text 
and  equation  (19)  above  by  substitution  of  the  various  quantities 
just  obtained. 
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It  has  been  explained  previously/  that  the  calculation  of  de¬ 
tonation  velocities  is  most  conveniently  carried  out  in  two  steps: 
an  ideal  detonation  velocity  is  first  computed  for  any  given  sub¬ 
stance,  then  the  actual  detonation  velocities  for  various  loading 
densities  are  calculated  from  the  ideal  value  with  the  assistance 
of  functions  previously  describedo 

By  combining  equations  (12)  and  (l4)  ve  obtain  for  the  ideal 


detonation  velocity: 


/  R  (j;  n 
V  M  / 


(3^) 


In  this  expression  nj|_  is  the  number  of  moles  of  gaseous  products 
formed  by  the  explosion  of  M  grams  of  the  explosive  material  and 
is  the  ratio  of  heat  capacities  averaged  over  the  gaseous 
products  at  constant  pressure  and  constant  volume  and  taken  at  the 
temperature,  '^21*  latter  is  to  be  calculated  from  equation  (13 

which  will  be  written  here  in  a  slightly  altered  form. 


_  _  Mfo -L'W'Wj  +T7 

Hf  denotes  heat  of  formation,  the  subscript  o  Indicates  the 


(35) 


original  explosive  and  the  subscript  j  refers  to  one  of  the 
products  of  the  explosion.  This  equation  has  been  obtained  from 
equation  (13)  by  introducing  the  substitlon: 
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and  the  approximation 

\  Cy-;  AT  r  (Qj  +  bj  'T, )  (37) 

"'t, 

vhlch  is  accurate  to  about  ifo  in  the  temperature  range  from  2000  to 

3000°K.  Numerical  values  of  the  coefficients  a.  and  b.  have 

J  ti 

been  given  in  Table  7.  The  sign  of  Ef  is  taken  positive  if  heat 
is  absorbed  in  the  formation  of  the  substance.  Any  quantity  of 
explosive  can  be  made  the  basis  of  calculation  provided  that  Hf^ 
and  the  various  nj’s  are  given  values  appropriate  for  that 
quantity.  All  refer  to  one  mole  of  the  respective  product. 

Since  it  Is  common  practice  to  tabulate  heats  of  formation  per  mole 
of  substance,  this  quantity  is  usually  the  most  convenient  basis 
for  calculation.  Tvo  small  departures  from  rigor  are  convenient 
vhen  using  equation  (35)  f03?  purposes  of  computation.  The  term 
involving  T2j_  vhlch  appears  in  the  numerator  makes  an  almost 
insignificant  contribution  and  its  value  for  T^^  equal  to  4000°K 
can  always  be  used  without  appreciable  error.  Likewise,  in  the 
denominator  the  term  involving  is  of  the  order  of  10^  of  its 
Companion  terms  so  that  an  error  of  10^  in  the  value  taken  for 
introduces  only  1^  error  in  the  value  for  T2i.  Since  the  total 
range  of  variation  of  values  of  for  substances  considered  In  this 
report  is  only  10^,  the  universal  use  of  the  mean  value  of  wil 

introduce  only  0,5^  uncertainty  in  the  value  for  Tgj^*  The  value 
for  P./2l'^'  Is  therefore  taken  always  to  be  0.80., 

The  calculation  can,  of  course,  be  carried  out  for  any  assumed 
distribution  of  products.  However,  only  one  set  of  products  is 
consistent  with  the  theory.  The  rules  to  be  followed  in  writing 
down  these  products  have  been  fully  discussed  in  the  text  and  will 
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merely  Toe  repeated  here.  These  are,  l)  only  solid  carbon,  CO, 

COg,  HgO,  Og,  Ng  (and  gaseous  sulfur,  SO2,  metals  or  metallic  oxides, 
insofar  as  the  corresponding  elements  are  present  in  the  explosix^e) 
are  formed,  and  2)  oxygen  reacts  quantitatively  to  form  CO,  any 
excess  then  reacts  quantitatively  to  form  H2O,  and  any  further 
excess  reacts  quantitatively  to  form  COg. 

Before  proceeding  to  Illustrate  the  method  of  computation  we 
must  discuss  one  further  convention  which  we  have  adopted  in  the 
use  of  equation  (3^)-  This  concerns  the  quantity  which  should 

be  ta.ken  at  the  temperature  calculated  from  equation  (35).  In  order 
to  obtain  the  value,  as  defined,  a  further  procedure  is  required 
which  Involves  the  use  of  tables  of  heat  capacities  at  various 
temperatures  from  which  appropriate  values  are  to  be  found  for  each 
of  the  product  substances,  after  which  a  summation  must  be  carried 
out  over  all  of  the  products.  This  procedure  is  fairly  tedious. 
However,  an  approximation  which  results  in  values  of  that  are 

uniformly  about  2%  too  low,  but  which  eventually  leads  to  values  of 
D  for  actual  conditions  from  which  this  error  has  been  eliminated 
is  very  simple.  Call  the  denominator  of  equation  (35)  d.  Then 

it  can  easily  be  shown  that 

—  _  0I+  £.8^; 

T:  -  d  +  ^58) 

As  stated,  the  use  of  this  approximation  in  equation  (3^)  results 
in  values  of  which  are  2^  low,  but  since  we  have  used  this 

procedure  consistently  for  all  calculations  the  error  is  compensated 
by  our  choice  of  covolumes  required  for  the  second  stage  of  the 
calculation  and  given  in  Table  6, 
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A  sample  calculation  is  nov  given  for  the  example. of  Nitroglycerin. 
The  calculation  is  based  on  one  mole.,  227  grams,  for  vhich  the 
proper  decomposition  equation  is: 


C3H5O9N3  =  3/2N2  +  5/2H2O  +  3CO2  +  1AO2.  Wi  =  7.25 

Hfo  =  -82,700  cal. /mole,  Hf(H20gg^g)  =  -57,870,  Hf(C02)  -  -94,240. 


n  -Hf  ■ 

J  J 

•T  a"’ 

3  3 

n 

1  1 

Wg  =  3/2 

8.31 

0.32  X  10"^ 

HgO  =  5/2 

-144,675 

19.65 

1.37  X  10“^ 

CO2  =  3 

-282,720 

30.48 

1.38  X  10“^ 

O2  =  1/4 

1.42 

0.05  X  10"^ 

-  N^.Hfi 

-  427,395 

NjAj  =59.86 

N,.b,=  3.12  X  10“^ 

j  J 

+  Hf^  = 

-  82,700 

344,700 

rn  „ 

344,700  +  300 

(59.86  +  3.‘ 

12  X  10“^  X  4  X  10-^) 

•^2i 

59.86  -  7 

.25  X  0.80 

+  3.12  X  10"^  Tg^ 

To.  = 

366,400 

53 .80  +  3.12  X  lO"-^  T 


Assume  T^.  =  5x10^;  Calculate  T2^  =  - -  =  “  5280°! 

AsauiDfcj  ^  ,  2i .  55,80  +  15.60  69.40 


Assume  =  5.20xlo5;  Calculate 


53.80  +  16.20 


366,400 
)  70.00 

5230°K  (accept.) 


d  =  70.00; 


70.00  +  2.8  X  7.25  89.6  _  n  on 

70.00  +  0.8  X  7.25  “  75.1  ~ 


Di  =  |n|0  ^ ^ . 5 . x.,1^2.Q  2  ^  ^q3  m/sec.  =  236O  m/sec. 


1)  From  Table  1 
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With  the  exception  of  nitroglycerin,  T.W.T.  and  the  salts 

of  metals,  the  calculations  for  all  other  substances  mentioned  in 

this  report  can  be  still  further  simplified  as  follows:  If  the 

formula  of  the  substance  be  written  as  C  H  0  then  q  ^  s  -  2q  +  r/' 

qrst  ' 

that  is,  there  is  at  least  enough  oxygen  to  convert  all  carbon  to  CO 
but  no  more  than  necessary  to  oxidize  all  hydrogen  and  carbon  to 
HgO  and  COg  respectively.  It  seems  almost  certain  that  the  majority 
of  the  compounds  likely  to  be  formed  from  these  four  elements  in  the 
search  for  new  and  powerful  explosives  will  answer  the  above 
description.  During  the  course  of  this  investigation,  it  was  found 
convenient  to  convert  equation  (35)  into  general  expressions 
involving  q,r,s  and  t.  These  are  somewhat  simpler  to  use  and 
furthermore  require  less  exercise  of  judgment,  or  rather,  attention 
to  rules,  and  do  away  with  the  necessity  for  reference  to  several 
of  the  tables.  In  the  interest  of  convenience,  their  use  is  strongly 
recommended  whenever  possible. 

When  r/2  for  a  substance  of  composition  CgH^OgW^ 

it  will  decompose  according  to  one  or  the  other  of  two  equations. 

A.  If  q  +  r/2  ^  s,  then  we  write: 

=  (t/2)W2  +  qCO  +  (s-q)H20'  +  (q-s+r/2)H2 

B.  If  q  +  r/2-^  s,  then  we  write: 

Vr°s^t  =  ^2°  +  (2q-s+r/2)C0  +  (s>q-r/2)C02 

Introducing  the  expressions  for  the  nj's  in  terms  of  q,  r,  s  and  t 
along  with  the  appropriate  numerical  quantities  taken  from  Table  7 
and  the  numerical  values  for  the  heats  of  formation  of  CO,  CO2  and 
H20j  then  proceeding  exactly  as  outlined  before  in  this  appendix  and 
as  Illustrated  for  nitroglycerin  one  obtains  the  two  equations: 
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^  Hfo  -  29,843q  +  l,222r,+  58,368s  +  9^8t 

^21=  ( 1 . 86q  f^06r+2 , 95  s+2 . 37t )  +  M  *  04q+0 .15^+0 . 25s+0  .lit)  .10"^ 


Q  +  p/2  £  s 


(ifO) 


IIIq  -  40,724q  -  4,2l8r  +  69,249s  +  948t 


"^21  (0.26q+1.26p+4.55s+2.37t)  +  ( -0.04q+0.15r+0. 253+0. lit ) .lO'^T^ 


For  "both  cases  the  equivalent  of  equation  (34)  becomes; 


Di  = 


+  1 


Vi. 

(t/2  +  r/2  +  q  ,  , 

— - - - .  8.31  f.  Tp.  .10  ^  ]  10^  m/sec. 

12q+r+l6s+l4t  j.  / 


The  use  of  one  of  these  equations  Is  Illustrated  by  the 
following  example,  a  calculation  for  cyclonlte. 

To  S3’’stematize  the  calculation  It  Is  convenient  to  consider 
both  equations  (39)  and  (40 )  to  be  of  the  form 


T  -  ^  ^1^1  .  _  P  y.  o  +• 

Xpi  _  ___ - - — ,  Xj^  -  q,r,s,u, 

^g^x^  +  10"^T2  2_hiXi 


Calculation  for 

cyclonlte. 

q  +  r/2  = 

=  3  =  s  =  3. 

Can  use  either 

(39)  or  (40) 

,  Arbitrarily  choose 

(40). 

^1 

hj^Xi 

q  =  3 

-122,172 

0.78 

-0.12 ■ 

r  =  6 

-  25,308 

7o6 

0.90 

s  =  6 

+415,494 

27.30 

1.50 

t  =  6 

5,688 

14.22 

0.66 

+421,182 

-147,480 

sura=49 . 86 

3.06 

-0.12 

sum  =  273,700 
+  Hfn  =  21,300 
295,000 


sum  =  2.94 


T  =  295,000 

2-  '49.86  +  2.94  X  10“3t, 
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Assume  Calculate  7^^=  Tif^^rxlrai  "  4680°ir 


Assu:ne  Tg^=K.65Klo5;  Calculate  12^=295^^  ^  MM 


4650^2 

acceot 


=  §1:^1±S-§.J.9  ^  1.2B 

^  63.5^  +  0.8  X  9 

D,  =  2^1^  /9  X  8.31  X  4.65  X  1.28^1^  «  2510  m/sec. 

^  1.28  222  / 


When  the  value  of  Dj_  has  been  obtained  for  a  given  substance  It 

becomes  a  relatively  simple  matter  to  predict  the  detonation  velocity 

for  various  loading  densities.  According  to  the  theory  given  In  the 

text,  D/Di  is  a  function  of  k  /M  vhere  D  Is  the  detonation  velocity/ 

at  the  loading  density,  ,  and  k  is  the  covolume  of  the  products 

resulting  from  decomposition  of  M  grams  of  the  explosive  at  the 

temperature  of  the  explosion.  We  begin  by  calculating  a  value  of  k 

at  300°.  The  sum  Is  taken  (over  all  products),  of  the  number  of 

moles  of  each  product  multiplied  by  its  covolume  as  given  In  Table  6. 

l/^ 

This  sum  is  then  multiplied  by  (300/T2i)  to  obtain  the  covolume 
at  the  temperature  of  the  explosion.  Further  multiplication  by 
the  desired  loading  density  and  division  by  M  then  gives  the  quantity 
vhlch  is  related  to  D/Dj_.  Numerical  values  of  D/D^  as  a  function  of 
numerical  values  of  k^t/M  are  given  in  Table  2  for  three  different 
values  of  as  defined  in  equation  (38). 

To  illustrate  we  consider  cyclonite  which  decomposes  by  the 
equation:  C^HgOgWg  =  3^2  +  3C0  +  3K2O  T2^  =  4650"  K  and 

=  2510  ra./sec.  Prom  Table  6  we  find  k^QQo  =  3(38.7)  +  3(38.7) 

+  3(25.8)  =  308.  Then  k^g^QO  =  308(300/4650)^/^  and  for  M  =  222 
and  ^,  =  1.20  we  obtain  k^/M  at  4650°K  =  O.665.  Observing  that 
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'fl  for  cyclonit e  is  1.26  and  consulting  Ta,ble  2  ve  find  in  the  ^rcl 
column  that  the  corresponding  value  of  D/Dj^  is  2.47.  Hence  ve 
predict  D  =  2.4?  x  2510  =  6l80  m/sec  at  f,  =*  1.20  g/cc.  If 

it  Is  desired  to  calculate  the  real  temperature  in  the  detonation 
vave,  the  knovledge  of  and  column  3  of  Table  5  make  this  very 
simple.  To  obtain  the  real  pressure  in  the  detonation  vave  one 
may  use  an  equation  vhlch  follovs  from  the  theory  developed  in 
this  report: 

P2  =  (x  +  x^e^)  (43) 

■  ^ 

After  calculating  k  as  discussed  previously,  one  may  take  the 
values  of  the  function  (x  +  x^e^)  from  the  last  column  of  Table  3 
and  thus  obtain  P  by  equation  (43)  in  vhich  R  =  82  cc.  Atm. 
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